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Abstract: The Gaussian expansion has been developed since early 80s as a powerful 
analytical method, which enables nonperturbative studies of various systems using 'pertur- 
bative' calculations. Recently the method has been used to suggest that 4d space-time is 
generated dynamically in a matrix model formulation of superstring theory. Here we clarify 
the nature of the method by applying it to exactly solvable one-matrix models with various 
kinds of potential including the ones unbounded from below and of the double-well type. 
We also formulate a prescription to include a linear term in the Gaussian action in a way 
consistent with the loop expansion, and test it in some concrete examples. We discuss a 
case where we obtain two distinct plateaus in the parameter space of the Gaussian action, 
corresponding to different large- solutions. This clarifies the situation encountered in the 
dynamical determination of the space-time dimensionality in the previous works. 
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1. Introduction 

The Gaussian expansion is a powerful nonpcrturbative method, which has been applied 
to condensed matter physics and statistical physics extensively. Even at the lowest order, 
which is sometimes referred to as the self-consistent harmonic approximation, it allows 
us to understand fundamental properties of various systems at least qualitatively. In fact 

there exists a systematic way to improve the approximation, which we refer to as the 
Gaussian expansion method, but it is also called in the literature under various names 
such as improved mean field approximation, delta expansion and so on. The method has 
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been worked out first in quantum mechanical systems [1, 2], where the expansion was shown 
to be convergent in some concrete examples [3], and it has been generalized to field theory 
later [4] . (The basic idea is also used in optimizing perturbation theory where the results 
depend on the renormalization scheme [2, 5, 6].) The most peculiar feature of this method, 
as emphasized by Stevenson [2], is that one obtains genuinely nonperturbative results — 
since never in the whole procedure does one attempt an expansion with respect to the 
coupling constant — and yet the required task is nothing more than familiar perturbative 
calculations based on Feynman diagrams. 

Recently this method has been applied to superstring/M theories using their matrix 
model formulations. In Ref. [7] Kabat and Lifschytz proposed to use the Gaussian approx- 
imation in the Matrix Theory [8], which is conjectured to be a nonperturbative definition 
of M-theory in the infinite momentum frame. Indeed Refs. [9] were able to reveal inter- 
esting blackhole thermodynamics from the dual strong-coupling gauge theory. An earlier 
application of the Gaussian approximation to random matrix models can be found in Ref. 
[10]. 

In Ref. [11] two of the present authors (J.N. and F.S.) applied the Gaussian expansion 
method to the IIB matrix model [12], which is conjectured to be a nonperturbative defini- 
tion of type IIB superstring theory in 10 dimensions. One of the most interesting questions 
in this model ^ concerns the possibility that the 4d space-time [17] appears dynamically 
accompanied with the SSB of the SO(IO) symmetry down to S0(4). From the path-integral 
point of view, this phenomenon may be caused by the phase of the fermion determinant 
[18], and Monte Carlo results support this mechanism [19] (Other possible mechanisms are 
discussed in Refs. [17, 20]). In Ref. [11] the same issue has been addressed analytically 
by using the Gaussian expansion method. The main idea was to consider various kinds of 
Gaussian action preserving only some subgroup of SO (10) and to identify the 'true vacuum' 
by comparing the corresponding free energy. Calculations up to the 3rd order showed that 
space-time preserving the S0(4) symmetry has the smallest free energy, and that the ratio 
of the extents in the four directions and the remaining six directions increases with the 
order. To our knowledge, this is the first time in history that the space-time dimensionality 
'4' is suggested from the nonperturbative dynamics of superstring theory. ^ 

A lot of effort has then been made to increase the order of the expansion. In Ref. [22] it 
was noticed that Schwinger-Dyson equations can be used to reduce considerably the number 
of Feynman diagrams to be evaluated. In Ref. [23] a computer code has been written in 
order to automatize the task of listing up and evaluating all the Feynman diagrams. With 
these technical developments, the order of the expansion has now been increased up to 
the 7th order, and the results strengthened the conclusion of Ref. [11]. The nature of the 
method itself has also been clarified. In Ref. [22] the Gaussian expansion method was 
interpreted as an improved Taylor expansion, and the importance of identifying a plateau 
in the space of free parameters in the Gaussian action was recognized. 

^The finiteness of the partition function has been proved in Refs. [13]. Various simplified versions were 
studied by Monte Carlo simulations [14, 15] and by the Gaussian expansion method [16]. 

^Recently another indication of this phenomenon has been obtained from the calculations of the 2-loop 
effective action around fuzzy-sphere solutions [21]. 
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Although these new results in matrix models obtained by the Gaussian expansion 
method are quite encouraging and deserve further investigations, we consider it equally 
important to know the method itself by applying it to a well-understood system and to 
try to improve or refine the method. In our previous work [24] we applied the method to 
the bosonic version of the IIB matrix model, where Monte Carlo results [14] are available, 
and the convergence of the method has been demonstrated. We also developed a new tech- 
nique to deal with the free parameters in the Gaussian action, which were conventionally 
determined by solving the 'self-consistency equations'. 

In this paper we attempt to gain more experience with the method by applying it to 
exactly solvable one-matrix models [25] , which have been studied intensively in the context 
of two-dimensional quantum gravity and non-critical string theory. The Gaussian expansion 
has been carried out maximally up to the 18th order and the results are compared with 
the known exact results. 

First we study the matrix model with various kinds of potential including the ones 
unbounded from below [25] and of the double- well type [26] . For the unbounded potential, 
it is known that a stable vacuum exists in the large-iV limit if the parameters of the 
potential satisfy a certain condition. In this case the Gaussian expansion converges and 
reproduces the exact results accurately. If the parameters are chosen such that there is 
no stable vacuum even in the large- limit, the Gaussian expansion docs not converge 
either. Thus the method captures correctly the critical phenomenon associated with the 
overflow of the eigenvalues. For the double-well potential, on the other hand, we find that 
the results are in reasonable agreement with the exact values at low orders, but they start 
to oscillate violently as we go to higher orders. 

Next we consider the 0^ matrix model. In this case the potential is always unbounded 
from below, but again a stable vacuum is known to exist in the large-A limit if the parame- 
ters of the potential are chosen appropriately. A new ingredient here is that the model does 
not have the Z2 symmetry unlike the (f)^ matrix model, and therefore we have to include 
a linear term in the Gaussian action. We formulate a prescription to deal with the linear 
term in a way consistent with the loop expansion. Since we have two free parameters in the 
Gaussian action in this case, identifying a plateau in the parameter space becomes more 
nontrivial than in the (j)'^ matrix model. Here we find the histogram technique proposed 
in our previous paper [24] to be quite useful, and we do obtain results converging to the 
exact values. 

Finally we reconsider the (f)^ matrix model with the double-well potential, in which the 
Gaussian expansion method in its simplest form seems to fail as mentioned above. Here we 
redo the analysis with the inclusion of a linear term in the Gaussian action. Two types of 
linear term are considered, and for each case we obtain a plateau with different free energy. 
When we use a SU(Ar) invariant but Z2 breaJcing linear term, we obtain a plateau which 
corresponds to putting all the eigenvalues into a single well. When we use a linear term 
which has the 'Z2 symmetry' but breaks the SU(Ar) symmetry down to SU(Ar/2)xSU(Ar/2), 
we obtain a plateau which corresponds to partitioning the eigenvalues equally into the two 
wells. The free energy for the latter is smaller in accord with exact results. This case 
therefore provides an example where we have more than one plateaus corresponding to 
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different large-A^ saddle-point solutions, but tlie true vacuum can still be determined by 
identifying the plateau which gives the smallest free energy. This supports the strategy 
taken to determine the dynamical space-time dimensionality in the aforementioned works 
in the IIB matrix model [11, 22, 23]. 

This paper is organized as follows. In Section 2 we study the 0^ matrix model with 
various types of potential. In Section 3 we study the (f)^ matrix model, where we discuss 
how to treat a linear term in the Gaussian action systematically. In Section 4 we revisit 
the matrix model with the double-well potential including a linear term in the Gaussian 
action. Section 5 is devoted to summary and discussions. In Appendix A we discuss general 
large-TV saddle-point solutions in the 0^ matrix model with the double-well potential. In 
Appendices B and C we present the details of the calculations we have done in Section 4. 

2. 0^ matrix model 

In this Section we apply the Gaussian expansion method to the (f)^ matrix model, which 
is exactly solvable and is known to exhibit nontrivial critical phenomena. We also give a 
brief review of the method in Section 2.2.1. 

2.1 Exact solutions 

The (j)^ matrix model we study in this Section is defined by the partition function 

Z = j d'''cj)e-^, (2.1) 

S = iVy tr 02 + 7v|tr (2.2) 

where cf) is a N x N hermitian matrix. The integration measure d^^ cf) is defined as d^'^cf) = 
Ua=i #"/\/2^, where are the coefficients in the expansion cp = (p'^T'^ with respect 
to the U(A^) generators normalized as tr (TT'') = {1/2)6"'^. By rescaling the variables 
as H^- \/\m\ (j), one finds that the system can be characterized by the sign of m and 
the effective coupling constant ^eflf = g/n^^- 

In order to discuss the property of this model, it is convenient to look at the eigenvalues 
of the matrix . Let us therefore diagonalize ^ as 

(f) = VAV^ , (2.3) 

where y is a unitary matrix and A = diag(Ai, • • • , A^r) is a real diagonal matrix. Integrating 
out the 'angular variables' V, one obtains the effective action for the eigenvalues Aj as 

N 

SeS = -J2^ri\Xi - \j\ , (2.4) 

Vix) = N (|x2 + , (2.5) 

where the second term in (2.4) comes from the measure d^'^cj). Thus the effective theory 
for the eigenvalues may be regarded as a system of N particles in the potential (2.5) with a 
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logarithmic repulsive force between each pair. As a fundamental object, we introduce the 
eigenvalue distribution 



(2.6) 



Let us also define the free energy F = —InZ and the observable (;^tr^^), which can be 
written in terms of p{x) as 



N 



-Mi 



dxx^p{x) . 



(2.7) 



(ii) two-lump 
phase 



(i)BIPZ phase 





c 

stable vacuum 



Figure 1: The phase diagram of the (p'^ matrix model. The 5 points labeled as A, B, • • •, E denote 
the set of parameters, for which we tested the Gaussian expansion method. 

In the large- limit the phase diagram of the model consists of three regions : (i) the 
BIPZ phase, (ii) the two-lump phase, and (iii) no stable vacuum, as depicted in Fig. 1. 

Let us first consider the case m > and > 0, in which the eigenvalue distribution 
p{x) is given by [25] 



p{x) = — ( -gx^ + ga^ + -m] a/ Aa^ — x"^ 
TT \2 2 J 



for |x| < 2a, where a is defined by 



m 



+ y^m? + l2g 



(2.8) 



(2.9) 



and p{x) = otherwise. Note that p{x) has a finite compact support [—2a, 2a] (Wc remind 
the reader that the large- A limit is already taken here). The absence of penetration into 
the region |x| > 2a is due to the fact that the potential (2.5) grows linearly as A goes 
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to infinity, while the distribution does not cohapse to the potential minimum due to the 
repulsive force between the N particles. The free energy and the observables are given by 



F{m,g) = 7V2 



i IniV - ^{ma' - l){ma' - 9) - ^ ln{2a') + O 0^ 



(2.10) 



\tr</>A = ia2(4- ma') , /i-tr </.A = — (3 - ma') . (2.11) 



-/V / 3 ' ' ' \N 

Let us set m to some positive value, and see what happens if one decreases g. Note that 
nothing singular happens at g = 0, and the above expressions remain valid for < as long 
as m' + 12^ inside the square root in eq. (2.9) is positive. Although the action becomes 
unbounded, the potential barrier grows with N fast enough to prevent the eigenvalues 
from overflowing. For g < — m'/12, however, the stable vacuum ceases to exist due to the 
overflow of the eigenvalues. At the critical point g = — m'/12, the non-analyticity of p{x) 
at the edge of the support {xq = ±2a) changes from \x — xqI^^' to \x — xqI^^'. This critical 
phenomenon plays a crucial role in the context of two-dimensional quantum gravity and 
non-critical string theory. 

Let us next set g to some positive value and see what happens if one decreases m. 
From (2.8) one finds that p"{0) changes from negative to positive at m = 2y^ , meaning 
that p{x) starts to develop a double-peak structure. This occurs before m reaches zero 
due to the repulsive force between the eigenvalues. From (2.8) one also finds that p(0) 
decreases as one decreases m and it finally becomes zero at m = —1^/g- Beyond this point 
the eigenvalue distribution p{x) will have two compact supports, and it is given explicitly 
by [26] 

p{x) = l.g\x\^{x^-Al)iAl-x^) (2.12) 

for < \x\ < A+, where A± is given by A^ = (|m| ±2^)/g and p{x) = otherwise. At 
m = —2y/g the distribution (2.12) reduces to (2.8). The free energy and the observables 
are given by [26] 

1 3 1, Im' 1, _/lniV\l 

,lniV----ln2--- + -ln, + 0(^jJ , (2.13) 

ltr</>A=^ + i. (2.14) 
N ^ I g^ g 

We call the phase that continues from the region m > and 51 > as the BIPZ phase, 
and the phase characterized by the two compact supports of the eigenvalue distribution as 
the two-lump phase. 

2.2 The Gaussian expansion method 

As we have seen in the previous Section, the simple one-matrix model exhibits nontrivial 
critical phenomena associated with the eigenvalue distribution. These phenomena are 
certainly nonperturbative in the sense that it cannot be seen by perturbative expansion 
with respect to g^s = g/m? . The aim of this Section is to apply the Gaussian expansion 
method to this model and to see if it reproduces these critical phenomena correctly. 



F{m,g) 






\m\ 








9 
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2.2.1 A brief review 

We will first illustrate the idea of the Gaussian expansion method using the cj)'^ matrix 
model as an example. Let us denote the quadratic term and the quartic term in the action 
(2.2) as ^2 and so that 8 = 82 + 84. Ordinary loop expansion can be formulated as 
follows : (i) consider the action 

8pt:{X) = j{82 + S4), (2.15) 

(ii) calculate free energy and observables as an expansion with respect to A up to some 
finite order, and (iii) set A = 1. Equivalence to the perturbative expansion with respect to 
9eS = g/rn^ can be easily seen by rescaling the variables as (f) i-^ \J A/|m| ^. The Gaussian 
expansion method simply amounts to replacing the action (2.15) by 



5GEM(t; ^) = X [{^0 + ^(-^2 - So)} + 54] , (2.16) 



A 

where is taken to be the SU(A'^) invariant Gaussian action ^ 

So = ^tr<^2 (2.17) 

with a real positive parameter t, which is left arbitrary at this point. Note that both (2.15) 
and (2.16) reduces to the original action (2.2) if one sets A = 1. The Gaussian expansion 
method can also be viewed as a loop expansion but with the 'classical action' (5*0 + 6*4) and 
the 'one-loop counterterms' [82 — 8q). Note, however, that the method cannot be viewed 
as an expansion with respect to any parameter in the original action unlike (2.15), except 
when one sets t = m, for which one simply retrieves the ordinary perturbation theory. In 
fact the freedom in changing the parameter t in the Gaussian action (2.17) plays a crucial 
role in this method. 

To make things more transparent, let us rescale 4> as (j) >—>■ so that the partition 

function takes the form 

Z = y"d^'0e-(^<=i+^'=* ) , (2.18) 
Sc\{t\ A) = 5o + A54 , 5c.t.(t; A) = \{82 - 80) . (2.19) 

Thus in the present example the Gaussian expansion can be viewed also as an expansion 
with respect to {8 — 80), the difference of the original action from the Gaussian action. 
However this is peculiar to the models with actions containing only quadratic and quartic 
terms, and it is not the case in general (See Section 3). In earlier formulation [7, 11] the 
expansion was first made with respect to {8 — So), and then it was reorganized in such a 
way that it becomes consistent with the loop expansion. Here we have chosen to take a 
shorter path. 

In actual calculations, the 'one-loop counter terms' can be incorporated easily by notic- 
ing the relation 

8ci{t;X) + 8c.t.{t;X) = 8ci{t + X{m-t);X) . (2.20) 
Gaussian action which breaks the SU(A'') symmetry is considered in Section 4.3. 
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Namely the calculation of a physical quantity in the Gaussian expansion proceeds in two 
steps; (i) obtain the A-expansion of the quantity using the 'classical action' Sci{t;X), (ii) 
shift the parameter t by A(m — t) and reorganize the expansion with respect to A. In 
general the first step can be done by ordinary Feynman diagrammatic calculations, where 
the use of Schwinger-Dyson equations reduces the number of diagrams considerably [22]. 
The large-A'^ limit can be taken, if one wishes, by simply retaining the planar diagrams 
only. (This results in considerable simplification in the matrix model applications, but 
it is definitely not the main point of the Gaussian expansion method itself.) In fact in 
the present case where we know the exact results for the 'classical action' Sci{t;X) [25], 
the first step can be done with little effort using softwares for symbolic manipulations like 
Mathematica. This allows us to test the Gaussian expansion method at very high orders 
(we went maximally up to the 18th order). 

Let us move on to the results obtained by the Gaussian expansion method. We compute 
the 'free energy density' defined by 

/= lim |^--ln7vl , (2.21) 

where the second term is subtracted in order to make the quantity finite in the large-iV 
limit. For illustrative purposes, we consider the case g = 1, m = (the point A in Fig. 
1), which was also studied in Refs. [22, 23]. Let us recall that the expansion parameter 
in ordinary perturbation theory is g^s = g/m? and the convergence radius is known to be 
1/12. Therefore the above choice of parameters corresponds to the strong coupling limit, 
which is of course beyond the reach of ordinary perturbation theory. 
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Figure 2: (Left) The free energy density / calculated by the Gaussian expansion method is plotted 

as a function of 1/t for g = 1, m = 0. Each curve corresponds to the order 3, 6, 9, 12, 15, 18. 
The horizontal line represents the exact value / = -3/8 - log(2/V3)/2 = -0.44692 • • •. (Right) 
The free energy density / extracted at the order 2 , • • •, 18 by the histogram technique from the 
results of the Gaussian expansion for various 1/t. (When we make a histogram, we used the region 
0.01 < 1/t < 1.) The 'error bars' represent the theoretical uncertainty estimated by the same 
technique. The order 1 result (/ = —0.423287) lies outside the displayed region. The horizontal 
line represents the exact result. 

In Fig. 2 (left) we plot the free energy obtained at the order 3, 6, 9, 12, 15, 18 as a 
function oi 1/t. As one can see also from the figure, the result of the Gaussian expansion 
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generically depends on the free parameter t in the Gaussian action (2.17). However, since 
t is a parameter which is introduced by hand, the result should not depend much on it in 
the region of t where the expansion becomes valid, if such a region exists at all. Indeed at 
sufficiently high orders we observe the formation of a plateau, meaning that there exists a 
certain range of t where the result becomes almost independent of t. Moreover the height 
of the plateau agrees very accurately with the exact result represented by the horizontal 
line. In old literature the free parameter (t in the present case) was determined in such a 
way that the result becomes most insensitive to the change of the parameter. This strategy 
was named as 'the principle of minimum sensitivity' [2]. The importance of the formation 
of a plateau has been first recognized in Ref. [22] . 

Given this new insight, we still have to develop some technique to identify a plateau and 
to extract its height in order to make a concrete prediction. Note also that the plateau is not 
completely flat but has some small fluctuations in practice. This causes some theoretical 
uncertainty in the prediction from the method, and it is desirable to be able to estimate 
its order of magnitude. In Ref. [24] we have developed a technique based on histograms. 
Fig. 2 (right) shows the results obtained for the present case. The 'error bars' represent 
an estimate for the theoretical uncertainty explained above. Glearly the result converges 
to the exact value within the first several orders. 



2.2.2 Unbounded potential (g < 0, m > 0) 

Let us move on to the g < 0, m > case, for which the potential becomes unbounded from 
below. As we discussed in Section 2.1, a stable vacuum continues to exist for g > — m^/12 
in the large- AT limit. In Fig. 3 (left) we show the free energy as a function of l/t at the 
critical point g = —1/12, m = 1 (the point B in Fig. 1). We observe the formation of a 
clear plateau. The results for g > — m^/12 were equally successful. It is noteworthy that 
the Gaussian expansion method can reproduce the large- solution even in the case of the 
unbounded potential. 
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Figure 3: (Left) The free energy density / is plotted as a function of l/t for 5 = — m = 1 
(critical point). Each curve corresponds to the order 3, 6, 9, 12, 15, 18. The horizontal line 
represents the exact value / — 7/24 — log 2 = —0.40148 • • •. (Right) The free energy density / is 
plotted as a function of l/t for 3 = — to = 1. Each curve corresponds to the order 3, 6, 9, 12, 
15, 18. 
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In Fig. 3 (right) wc show the results for g = —j^, m = 1 (the point C in Fig. 1), which 
is shghtly below the critical point. At each order we see a plateau-like region, but the slope 
in that region increases with the order (Note the striking difference from the figure on the 
left). Thus the Gaussian expansion method seems to know the absence of a stable vacuum 
with this choice of parameters. This example also confirms the importance of the plateau 
formation in the method. 

It is interesting to sec how the plateau formation ceases to occur when we cross the 
critical point. This is reflected in the peculiar situation at the critical point. The plateaus 
in Fig. 3 (left) do not have the tiny oscillations which become visible upon magnification 
in Fig. 2 (left). In fact in Fig. 3 (left) there is only one stationary point at each order we 
studied, and the free energy at that point coincides with the exact value at the orders > 2. 
The stationary point disappears as soon as we cross the critical point. 



2.2.3 Double-well potential {g > 0, m < 0) 

Let us then consider the case with the double- well potential. The system undergoes a phase 
transition from the BIPZ phase to the two-lump phase at m = — 2y^. Fig. 4 (left) shows 
the free energy as a function of 1/t for g = 1, m = —1 (the point D in Fig. 1). Although the 
chosen parameters are still well within the BIPZ phase, we observe an oscillating behavior 
which becomes more and more violent as the order is increased beyond the order 9. This is 
also reflected in Fig. 4 (right), which shows the results extracted by the histogram technique 
at each order. The oscillation becomes even worse as we go to further negative m. 
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Figure 4: (Left) The free energy density / is plotted as a function of l/t for g = 1, to = — 1. 
Each curve corresponds to the order 3, 6, 9, 12, 15, 18. The horizontal line represents the exact 
value. (Right) The free energy density / obtained by the histogram prescription at the order 2, 
18 for = 1, TO = —1. (When we make a histogram, we used the region 0.01 < l/t < 1.) The 
order 1 result (/ = —0.846574) lies outside the displayed region. The horizontal line represents the 
exact result / = -0.933705 • • •. 
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Thus we find that the Gaussian expansion method in its simplest form does not work 
in the case of the double-well potential. (A similar result was obtained in the quantum 
mechanics of a particle in a double- well potential [1].) This is understandable if we recall 
that the method amounts to expanding the theory around ^ = 0, which is a local maximum 
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of the double-well potential. We will reconsider the present example in Section 4 including 
a linear term in the Gaussian action. 



3. 0^ matrix model 

In this Section we will test the Gaussian expansion method in the (f)^ matrix model, which 
is also exactly solvable. In this model the potential is always unbounded from below, but a 
stable vacuum exists in the large- limit when the parameters in the potential are chosen 
appropriately. The situation is thus similar to the 0^ model with the unbounded potential, 
but the main issue here is that we have to include a linear term in the Gaussian action 
due to the absence of the Z2 symmetry. We will treat it in a way consistent with the loop 
expansion, and check the result against the exact solution. As a result of having a linear 
term, we have two free parameters, which makes the identification of a plateau somewhat 
more nontrivial than in the previous Section. The histogram technique turns out to be 
quite useful. 

3.1 Exact solutions 

The (f)^ matrix model is defined by the partition function 



Z 



j d'^'ct^e-^, (3.1) 



S = ivitr02 + iv|tr(/)^ , (3.2) 

where (pis a. N xN hermitian matrix, and the measure d^^ 4> is the same as in the matrix 
model. Although the action (3.2) is unbounded from below, the theory is well-defined in 
the large-iV limit for cp- < g\ = 

The exact results for the free energy and observables are given by [25] 



F{g) = N- 



1 ll + 2a 1(7(30-2-^6(7-^2) ^ flnN 
m A/ -|- - m — — — ^ . „ -|- O 



2 2 2 3(l + (7)(l + 2(7)2 V N 



1 , A 1 + 3cr / 1 , ,2\ 1 + 3C7 



, (3.3) 

where a is the largest solution of the equation a{l + a){l + 2a) = —2g^. The above 
expressions are valid for 1^1 < Qc- For \g\ > gc, on the other hand, the system does not 
have a stable vacuum even in the large- iV limit. In what follows we assume g >0 without 
loss of generality due to the duality under the sign flip of g. 

3.2 Systematic treatment of a linear term in the Gaussian action 

Let us denote the quadratic and cubic terms in the action (3.2) as S2 and Ss so that 
the action reads 8 = 82 + S3. Similarly to (2.16), the Gaussian expansion amounts to 
considering the action 

SGEuit, A) = ^ [[So + \{S2 -So)]+ Ss] . (3.5) 
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Since the action (3.2) does not have the Z2 symmetry unUke the 0^ matrix model, the Sq 
in (3.5) should include a linear term as 



So = Sg + Si,, 

Sg = ^tticfy^, Sl = Nhtr(l) , (3.6) 

where the real parameters t and h are arbitrary at this point. Free energy and observables 
are calculated as an expansion with respect to A, and A is set to one eventually. This 
procedure yields a loop expansion with the 'classical action' {Sq + S3) and the 'one-loop 
counter term' (^2 — Sq)- 

As in the (j)'^ matrix model, let us rescale (f> as cj) >—>■ so that the partition function 
takes the form 

Z = j d^>e-(^'='+^'= *-) , (3.7) 
Sa{t, h-X) = Sg + ^Si^ + VXSs , S,.t. {t, h; A) = \{S2 - Sg) - VA-Sl ■ (3.8) 

Thus the Gaussian expansion in the present case cannot be viewed as an expansion with 
respect to {S — So) unlike the 0^ matrix model. The 'counter terms' ^ct. in (3.8) can be 
easily implemented by using the relation 

5ci(t, h; A) + 5c.t.(i, h; A) = Sd{t + A(l -t),h- Xh; X) . (3.9) 

Namely the calculation of a physical quantity in the Gaussian expansion proceeds in two 
steps; (i) obtain the A-expansion of the quantity using the 'classical action' Sc\{t,h; X), 
(ii) shift the parameters t and h by A(l — t) and {—Xh) respectively and reorganize the 
expansion with respect to A. 

As usual we eliminate the linear term in Sci{t,h;X) by shifting the variables as = 
^ + ^1, where the parameter a should satisfy 

ga^ + ta + h = . (3.10) 

The 'classical action' then becomes 

Sait,h;X) = C + Nl-i tip + VXN\gti , (3.11) 

where C and i are given by 

C = N^j Qa^t + ha + ^ga^^ , (3.12) 
i=t + 2ga. (3.13) 

Thus the problem reduces to the conventional pcrturbativc expansion in the (f)^ matrix 
model. Note that the constant term C in (3.11) contributes to the 'order-(— 1) term' in the 
A-expansion of free energy. 
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14000 




-0.3895 -0.389 -0.3885 

/ 

Figure 5: (Left) Free energy density at the order 5 is plotted as a function of 1/t and h. (Right) 

The histogram of the free energy density. (When we make a histogram, we nscd the region 0.02 < 
1/i < 2, — 3 < ft < 5. The bin size is taken to be 10~^.) The vertical line represents the exact 
result / = -0.389176 • • •. 




Figure 6: The same as Fig. 5 but for the order 15. 



Eq. (3.10) has two solutions a = ^(— i ± \Jt'^ — 4^^), which correspond to t = 

iby^ — 4g'/i respectively. Since t is the coefficient of the quadratic term in cq. (3.11), 

it should be positive in order for the Gaussian expansion to be well-defined. This forces us 

to take ^ 

a = — (-t + V't^-45/i,). (3.14) 
^9 

Considering that 1/t has the meaning as the tree-level propagator, we shall plot the result 
of the Gaussian expansion obtained for various t and ft. as a function of 1/t and ft. For 
each t and ft satisfying P + \gh > 0, we obtain results corresponding to t > and t < 0. 
It turned out, however, that the plateau mainly develops in the branch corresponding to 
t > 0. We therefore plot the results for t > in what follows. 

Here we show our results at the critical point = g^ = In Fig. 5 (left) we plot 

the free energy density defined by (2.21) as a function of l/t and h at the order 5. We 
observe a clear plateau. On the right we show the corresponding histogram, which has a 
sharp peak at the exact value. Fig. 6 shows the results at the order 15. The plateau extends 
and as a result the peak of the histogram becomes higher. (In the present case we were not 
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able to obtain a reasonable estimate for the uncertainty because the plateau turned out to 
be too flat. This seems to be a peculiar feature at the critical point.) We have checked that 
observables such as (^^tr^) and {^trcj)'^) also agree with the exact results very accurately. 
These results support the correctness of the way we treated the linear term in the Gaussian 
action. 

4. (j)'*' matrix model with the double-well potential, revisited 

In this Section we will revisit the matrix model with the double-well potential, where 
the simplest version of the Gaussian expansion seems to fail (See Section 2.2.3). Here 
we would like to perform the expansion including a linear term as we formulated in the 
previous Section. In fact we are going to consider two kinds of linear term. One is the 
same as in the cf)^ matrix model, and the other is a linear term with a 'twist' which breaks 
the SU(A/") symmetry down to S'U{N/2)xSU{N/2). For each case we obtain a plateau 
corresponding to different large-A?^ solutions. The first one corresponds to putting all the 
eigenvalues into one well, and the second one corresponds to partitioning the eigenvalues 
equally into the two wells. In accord with the exact results, the free energy is smaller for the 
latter, which corresponds to the 'true vacuum'. The present example supports the validity 
of the strategy adopted in the dynamical determination of the space-time dimensionality 
in the IIB matrix model [11, 22, 23]. 

4.1 Single-lump solution in the double- well potential 

In fact the 0^ matrix model with the double-well potential has a one-parameter family 
of large- AT solutions corresponding to different partitioning of the N eigenvalues into the 
two wells (See Appendix A). The solution (2.12) corresponds to the equal partitioning, 
and it has the smallest free energy. For asymmetric partitioning, the free energy becomes 
larger due to the repulsive force and due also to the increase of the potential energy. These 
asymmetric largc-A^ sohitions are stable, however, against tunneling of the eigenvalues 
because the potential barrier grows linearly with N. The situation is analogous to the 
existence of a stable large-A^" solution in spite of the unbounded potential discussed in 
Section 2.1. 

Here we will consider the extreme case where all the eigenvalues are put into one well. 
The exact results can be readily obtained by applying the technique used in Ref. [25]. The 
eigenvalue distribution has the form 



for < X < Bj^ and p{x) = otherwise. The edges of the support are given by 




(4.1) 




(4.2) 



and the parameter a is defined as 



cr = 



3|m| + 2^m^ - 15^ 
5 



(4.3) 
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The free energy and observables are obtained as 



1 3 1 So- 
- In iV h - In 

2 8 2 4 

(5 2 m?\ \m\ 
3 ~ 45^7 5a 



43 



+ 



180 5 
InA^ 



N 



-tr( 



g 5a — \m\ ' \ A?" 



1 



-tr( 



16 



41 \m\ 
45T^^ 



1 mr 
15 7" 



- 1 



5cr — \m\ 



(4.4) 
.(4.5) 



As can be seen from (4.3), the single-lump solution exists for < g/m^ < 1/15, while the 
Z2-symmetric solution (2.12) corresponding to the 'true vacuum' exists for < g/m^ < 1/4. 
It is understandable that the critical g/m? is larger for the latter (the potential wells become 
shallower as g/m? increases). 

4.2 Gaussian expansion method with a linear term 

We will see that the single-lump solution (4.1) can be reproduced by the Gaussian expansion 
method if we include a linear term. The formulation is quite analogous to the cf)^ matrix 
model, so the description will be brief. Instead of (3.5) we consider the action 



SoEuit, h;\) = i -t- A(S'2 - So)! + S4 



(4.6) 



where includes the linear term Si, as in (3.6), while S2 and S4 are the same as in Section 
2.2.1. Rescaling the variables as cf) V\(t>-, we obtain 



S,.tXt, h; A) = X{S2 - Sg) - VxSi 
The 'one-loop counter terms' can be incorporated by using the relation 
Sci{t, h; A) -|- S'c.t.(ij h; A) = Sc\{t + A(m — t),h — Xh; A) . 



(4.7) 
(4.8) 

(4.9) 



After eliminating the linear term in Sci{t, h; A) by shifting the variable as (f) ■ 
'classical action' becomes 

Sci{t,h;X) = C + N^i ti (f>^ + ^/XN^gtI(f)^ + XN^gtT ^'^ , 



<+^l,the 



(4.10) 



where C, i and g are given by 



C = N- 



1 a 
A4 



(3/1 + ta), 



t = t + 3ga , 
9 = ^ga . 



(4.11) 

(4.12) 
(4.13) 



Thus the problem reduces to the conventional perturbative expansion in the one-matrix 
model with both c/)^ and cp'^ interactions. As in the previous cases, we can utilize the large- A^ 
solution, which can be obtained by the same method [25]. See Appendix B for the details. 
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The shift parameter a should satisfy 

go? + ta + h = Q, (4.14) 
which has three solutions. For f > they arc given by 

18V3^ ^+ e ^- J ' 



02 



as 



(e-W3GV3 _ ,iV3Gry3) ^ (4.15) 



while for t < they are given by 



1 ../3(cV3_gy3-) ^ 



where we have introduced G± = ±9-^h+ \/^lgh?' + 12fi. Among the three branches, two 
of them (a2 and 0:3 for t > Q, a\ and a2 for t < 0) are mapped to each other under the 
sign flip of h, while the remaining one is symmetric. This structure is a consequence of the 
Z2 symmetry of the model. Requiring i to be real positive we are left with the following 
four cases. 

a = ai for t > , 



a = ai for t < 0, h < - ^ ^\t\^/'^ , 

a = a2 for t < 0, h> ^f^W'^ , 

a = ai,a2 for t< 0, -J^^IA^^^ <h< J^W"^ ■ 



(4.17) 



As in the (jy^ model, we plot the result of the Gaussian expansion as a function of 1/t and h. 
In what follows we show only the branch which includes the plateau in the /i > regime. 

We present the result for the free energy at 5 = 1/15 and m = — 1 (the point E in Fig. 
1), which is the critical point for the existence of the single-lump solution. In Fig. 7 (left) 
the free energy density / at the order 3 is plotted as a function of 1/t and h. The figure 
on the right shows the corresponding histogram. (Here again we were not able to obtain 
a reasonable estimate for the uncertainty because the plateau turned out to be too flat.) 
The formation of a plateau is clearly observed and the histogram has a sharp peak at the 
exact result. Fig. 8 shows the results at the order 6. The plateau extends to larger area, 
and consequently the peak of the histogram becomes higher. 

4.3 Gaussian expansion method with a twisted linear term 

In this Section we attempt to reproduce the 'true vacuum', which corresponds to parti- 
tioning the eigenvalues equally into the two wells. The failure of the Gaussian expansion 



-16- 







-3.8055 



3.8045 



Figure 7: (Left) Free energy density at the order 3 is plotted as a function of 1/f and h for 

g = 1/15 and m = — 1. (Right) The histogram of the free energy density. (When we make a 
histogram, we used the region 0.02 < 1/i < 2, —5 < h < 5. The bin size is taken to be 10~^.) The 
vertical line represents the exact result / = —3.80491 • ■ •. 



8000 




Figure 8: The same as Fig. 7 but for the order 6. 



method in Section 2.2.3 may be attributed to the fact that the A-expansion with the action 
(2.16) corresponds to expanding around the classical solution (p = 0, where the double-well 
potential becomes locally maximum. (When we say 'classical', we regard the parameter A 
as the Planck constant h.) Here we will take the Gaussian term 5*0 in the action (2.16) in 
such a way that the classical solution becomes proportional to 

J = (73 ® lLjv/2 = diag(l,-y,l, -l,---,-l ) , (4.18) 

N/2 N/2 

where and henceforth aa denote Pauli matrices. This can be achieved by considering the 
action (4.6) with a 'twisted linear term' given by 

Sl = Nhtr{J(l)) . (4.19) 

This term breaks the SU(A^) symmetry down to SU(A^/2)xSU(iV/2), but the total action 
(4.6) has the 'Z2 symmetry' 



(f) —T(f)T^ , T = c7i (8) I7V/2 



(4.20) 



We expect that other large-A^ solutions described in Appendix A can be reproduced by 
considering the asymmetric twist 

/ = diag(l,---,l, -l,---,-l ) , nL + nR = iV (4.21) 

instead of (4.18). 

After eliminating the linear term by shifting the variable as (f) = (f)+^J, the 'classical 
action' becomes 

Sci{t, h;\) = C + N^i ti 4>^ + N^g tr ( J0J(/>) + VXNg tr ( J^^) + XN^g tr , (4.22) 
where C, t and g are given by 

C 
i 
9 

By decomposing the matrix ^ as 

4 

4> = ^<ya®Va , (4.26) 
a=l 

where (T4 is a 2 x 2 unit matrix and (pa are Ar/2 x Ar/2 hermitian matrices, the quadratic 
term in (4.22) becomes 

5quad = N{i-ga^)iT{{ipif + {^2?] 

+N{i + g a^) tr {{yp^f + (.^a?} ■ (4.27) 

In order to make the Gaussian expansion well-defined, we therefore have to require t—g o? > 
0, which means 

t + ga^>0. (4.28) 

The equation for a is the same as (4.14), so the solutions are (4.15), (4.16). Combining 
(4.28) and (4.14), one finds that a and h should have opposite signs. Thus the allowed real 
solution for a is determined uniquely in the present case as 

a = ai for t > 

a = ai for t < 0, h < (4.29) 
a = a2 for t < 0, h > . 

Note that t and h can take any real values except for t < and h = 0. As before we plot 
the results obtained for various t and h as a function of — gcP') and h. In the present 
case we have only one branch. Since the figure is invariant under the reflection h ^ —h 
due to the Z2 symmetry, we will restrict ourselves to /i > in what follows. 



t + 2g o?, 
ga . 



(4.23) 

(4.24) 
(4.25) 
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Figure 9: (Left) A contour plot of the free energy density as a function — ga^) and h at 
the order 1 for m = — 1, g = 40 contours for / = —5.5 + 0.02 x n (n = 1, • • • , 40) are plotted. 
The blank region near the right-up corner means that the calculated free energy density does not 
fall in the specified region of /. (Right) The same as in the figure on the left but for the order 3. 
The plateau has grown diagonally from the left-top to the right-bottom. The meaning of the blank 
regions is the same as in the left figure. 



quantity 


/ 






order 1 


-4.82843 


15.8906 


241.361 


order 3 


-4.82861 


15.8899 


241.346 


exact 


-5.14859 


15 


240 



Since the one-matrix model with the SU(A^) breaking terms (4.22) cannot be solved 
exactly with the known techniques, we have to perform diagrammatic calculations explic- 
itly unlike the previous cases (See Appendix C for the details). We have performed the 
calculations up to the order 3. In Fig. 9 we show contour plots of the free energy density 
as a function of l/(t — ga^) and h at the order 1 and 3form = — 1,^ = ^ (the point E 
in Fig. 1; the same as in Section 4.2). 

Since our analysis here is restricted 
to low orders, where the use of the his- 
togram method is not highly motivated, 
we evaluated the free energy at the ex- 
trema on the plateau as in Refs. [11, 22, 
23]. Although there are more than one 
extrema at the order 3, the difference in 
the results is negligible within the preci- 
sion we are discussing. The observables 
(;^tr^^) and (^^tr^^) are also evaluated at the extrema of the free energy as in Refs. 
[11, 22, 23]. The results arc shown in Table 1. We see a reasonable agreement with the ex- 
act results already at the order 1. The trend of improvement seen at the order 3, however, 
turns out to be tiny. This might be related to the fact that the twisted linear term breaks 
the full SU(A^) symmetry, which is expected to be restored at higher orders. 

Let us emphasize that the results in Sections 4.2 and 4.3 are obtained for the same 
action with the same parameters. Depending on the choice of the Gaussian action, we have 
seen plateaus with different free energy. Each plateau corresponds to a large- saddle-point 
solution, and the true vacuum can be identified with the plateau which has the smallest 
free energy. Thus our results support the strategy adopted in the dynamical determination 
of the space-time dimensionality in the IIB matrix model [11, 22, 23]. 



Table 1: The results of the free energy and the ob- 
servables at the order 1 and 3. 
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5. Summciry and discussions 



In this paper we have tested the Gaussian expansion method in the one-matrix model, 
which is exactly solvable and is well-understood in the context of 2d quantum gravity and 
non-critical string theory. 

One of the most important findings is that the method is able to reproduce large-A^ 
solutions which are not stable at finite N, but stabilize only in the large- A'^ limit. Examples 
we have seen are the vacuum in the unbounded potential (as in the 0^ model with g < 0, 
m > and the 0^ model) and the single-lump solution in the double-well potential (0^ 
model with g > 0, m < 0). This enables us to understand the situation encountered in the 
study of the TIB matrix model, where we obtained many plateaus with different space-time 
dimensionality. We may naturally expect that each plateau corresponds to some large- 
saddle-point solution and that the true vacuum can be identified with the plateau which 
gives the smallest free energy as in the case of the double-well potential. In fact this was 
assumed implicitly in the dynamical determination of the space-time dimensionality, but 
now we have provided a concrete example where the statement holds. 

On the technical side, we would like to emphasize the importance of the plateau for- 
mation in this method. We have observed in the case of the unbounded potential that the 
plateau formation is quite sensitive to the stability of the vacuum. We have also formulated 
a prescription to include a linear term, which was shown to work in various examples. In 
these cases we have to deal with two free parameters, which makes the identification of a 
plateau more nontrivial. Here we find the histogram technique to be quite useful. 

To conclude, we hope our results clarified the nature of the Gaussian expansion method 
and confirmed its usefulness in particular in matrix model applications. We expect that 
the method can still be refined or extended in many different ways, thus allowing us to 
study various systems which are not easily accessible by other means. 
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A. General solutions in the double-well potential 

In this Appendix we discuss the general large-A^ solutions in the double-well potential 
(m < 0, y > case of the cp'^ matrix model). Each solution corresponds to different 
partitioning of the N eigenvalues into the two wells, say, nj^ in the left well and tir, in 
the right well, where ni, + = N. In the large- A?^ limit, the support of the eigenvalue 
distribution is given by L = [—B', —A'] U [A, B]. As a parameter which characterizes each 
solution, we introduce 




(A.l) 
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which takes the values = 0, 1, oo corresponding to (riL, riR) = (0, N), {N/2, N/2), {N, 0), 
respectively. Since the solution for l/v is Z2-equivalent to that for v, we may restrict 
ourselves to < 1. Introducing the variables 



(A.2) 



the saddle-point condition reads 

(s — i's'){\m\ — — s'^ + ss') + '—^—ss{s — s) 



0, 



1 - V 



J\2 



+ 



2 



(|m| -s'- + ss') 



(|m| 



s" + ss'){s + s'){{2 - u)s + {2u - l)s'} = g{l + , (A.3) 



which determines s and s' for any u, and the eigenvalue distribution has the form 



p{x) 



J 2^ + - s')} ^{ x ^ A){B - x){x + A')ix + W) for x G [A, B] 
1 2¥ [-9^ - Vais - s')] ^{A -x){B- x){-x - A'){x + B') for x € [-B', -A'] . 

(A.4) 

The equations (A.3) are not solvable analytically for general v except for 1/ = 0, 1, at which 
one obtains the analytic solutions (2.12) and (4.1) respectively. 



B. Exact solution for the model with both 0^ and 0^ interactions 

In Section 4.2 the problem reduced to the A-expansion of the free energy defined by 



exp[-iV 



exp 



-ivihr(/)2 - ^/^N^gii4)^ - XN^gti(t>^ 



(B.l) 



The one-matrix model with both 0^ and interactions can be solved in the large- AT limit 

using the method of Ref. [25]. Assuming that the eigenvalue distribution p{x) defined by 
(2.6) has a finite support [2a, 26] in the large- A'^ limit, it takes the form 



p{x) = ^{Xgx^ -Ax- B)^/{x-2a){2b-x) 
27r 

A = -\g{a + h)- VXg , 



B = -i 



-Xg{b - a)^ - \g{a + b)^ - V\g{a + b) . 



In order to simplify the expressions, we introduce 

g = pK, g = i'2k, a + b 

where the parameters s and A are determined by 
-2s 



ks , {b — a 



|2 = iA 
t 



A 



(1 -I- Xk'^S + X^HiK^S^) , 



1 -|- 2>Xks 

-2 
1 + 12Ak 
1 



I + -X{k + 2k^ + 12Xk^)s^ + X^{bK + k^)k^s^ 



4 2 4 



(B.2) 

(B.3) 
(B.4) 



(B.5) 



(B.6) 
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In terms of s and A, the free energy is given as 

jr=LnAr-^ + iln^ + -X^h' + -A^^S^ + h^KK^'s^ + -(1 + 2XKh + 3X\kh')A 
242A2 3 4 16 

5 1 

+— AkA^ + —Xk^s{1 + Xk^s + A\k^s^)(1 + 3Aks)A^ 

+— A f K + -^2(2 + 9Aks) + AK/t2s(4 + 15Aks) ) A^ + ^-A^^^A'' 
256 V 3 ^ ' ^ 7 2048 

Although we cannot solve (B.6) with respect to s in a closed form, we can obtain s order 
by order in A by iteration, which is sufficient for our purpose. Up to 0(A®), the parameter 
s and the free energy T are given respectively as 

s = -2 + 6A(3k - 2k^) - 4A^(45k2 - mnk^ + 32k'') 
+10A^(189k^ - 786^2 + 712kk^ - 168k^) 
-12A^(1701ac'' - mh2K^k^ + 2mUK^k'^ - lUQAKk^ + 2048^^) 
+28A^(8019k^ - 85482k^k2 ^ 228240k3ac^ - 230032k2k6 + 95488kk8 - 13728k^°) 
-8A^(312741k^ - 4712580k^k2 + 18297792k''k^ - 28301808ac^ac^ + 20097664/«2k8 
-6537600kk^° + 786432k^2) + 0{X') , (B.8) 

J^=JlniV + ^ln^ + ^A(3k - Ak^) + X^ ( -^k^ + Qkk^ - 
2 Z 2 \ " 3 

+A3 [^^3 _ 54^2-2 ^ g8^~4 _ 

+A4 f-i^K^ ^ 504^3-2 _ 1236^^2 /c^ + 856^^6 _ 512 
Y 8 3 

(729 91 52 
- 4860k^k2 + 19584^3-4 _ 24608^2 + 11680kk^ k^^ 
5 5 

+A6 / _5^^6 ^ 48114^5-2 _ 288468ac''k^ + 567816k^k^ - 469472^2^8 



+169152^/^1° - ^^12^ + 0{X') . (B.9) 

C. Diagrammatic calculations with a twisted linear term 

In this Appendix we describe the diagrammatic calculations in the presence of the twisted 
linear term. As we discussed in Section 4.3, the problem reduces to the A-expansion of the 
free energy r) defined by 



-N^H tr (j)^ - N^r tr {J(j)J4>) 
-VXNgtv{J^^) - XN^gtvct)"^ 



(C.l) 
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where the parameters ^ and r are related to those in (4.22) as ji = i and r = go? . The 
Feynman rules are given in Fig. 10. Instead of evaluating all the diagrams contributing 
to J^(//, r), we use the Schwinger-Dyson equations to reduce the number of diagrams to 
be evaluated. Here we extend the idea of Ref. [22] to the case including a linear term. In 
what follows we restrict ourselves to the large- A?^ limit, so that we have to consider planar 
diagrams only, but the method itself is applicable to finite N as well. 
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Figure 10: The Feynman rules for the A-expansion of (C.l). 



C.l Full propagator, tadpole and 2PI vacuum diagrams 

As the fundamental quantities in the present calculation, we consider the full propagator 
(the connected two-point function) and the tadpole (the one-point function), which can be 
written respectively as 

{<t>ij(t>kl)c = {<f>ij(f>kl) - {(l>ij){<f>kl) = j^{c + dJiiJjj)5ii5jk , 

(cPij) = VXTJiAj (C.2) 

due to the SU(iV/2)xSU(Ar/2) symmetry. At the leading order, c, d and T are given by 

c=^^ + 0{X) , d=^^ + 0{X) , T = -g ^ ^^~^ +0(A). (C.3) 

In what follows we will use the Schwinger-Dyson equations for the full propagator and 
the tadpole to reduce the calculation of the free energy to that of two-particle-irreducible 
(2PI) planar vacuum diagrams. By 'two-particle-irrcducible' we mean that the diagram 
cannot be separated into two parts by removing two propagators. For later convenience 
we use the full propagator for internal lines, and denote the sum of those diagrams as 
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N'^G{c,d). For example N'^G{c,d) consists of 5 diagrams shown in Fig. 11 up to the 2nd 
order, and the results of the first two diagrams are given by ^Xg^N'^{c^ + (P + 2(?d + 
2cd^) and —^XgN^c^ respectively. The important point is that there are much less two- 
particle-irreducible (2PI) vacuum diagrams than general vacuum diagrams that need to be 
considered in order to obtain the free energy J^{fi, r) directly. 



N G(c,d) 








+ 0( X ) 

Figure 11: The 2PI planar vacuum diagrams up to O(A^). The internal line with the rectangular 
blob stands for the full propagator. The two diagrams in the first line are of 0(A), while the three 
diagrams in the second line are of O(A^). 

C.2 Derivation of the Schwinger-Dyson equations 

Here we will derive the Schwinger-Dyson equations, which allow us to obtain the full 
propagator and the tadpole order by order in A. 



NBj. 



Figure 12: The full propagator can be expressed as a geometric scries. The round blob represents 
the radiative correction to the quadratic term, which is denoted by NBij as indicated in the second 
line. 

Let us note first that the full propagator can be expressed as a geometric series as 
depicted in Fig. 12. The round blob, which we denote as NBij, stands for the radiative 
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Figure 13: The first two lines show the diagrams which contribute to the parts of Bij and T 
which depend only on T but not on c or d. The round blob attached at the ends of the double-line 
represents the tadpole TJaSij as stated in the last line. 




Figure 14: The full result for Bij (c, d, T) can be obtained from Bij (c, d,T = 0) by dressing the 
3-point vertices with the tadpole. 

correction to the quadratic term. One can solve this relation for B^j as 

It is convenient to decompose B^j as 

B,j = Bij{T) + Bij{c,d,T) , (C.5) 

where Bij{T) contains only tadpoles and hence do not depend on c or d. One can easily 
see that Bij{T) consists of 5 diagrams shown in Fig. 13, which yield 

Bij{T) = {-2X~gT - X^gT^)(2 + JuJjj) • (C.6) 

The second term in (C.5) can be obtained from the sum of 2PI diagrams G{c,d) as fol- 
lows. Note first that Bij(c,d,T = 0) can be obtained from G{c,d) by removing one full 
propagator, i.e., 

Bij{c,d,T = 0) = 2—G{c,d) + 2—Gic,d)JiiJjj ■ (C.7) 
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N G(c,d) 



T(c,d,T= 0) Jji 




Figure 15: Except for the lowest order diagram shown in the upper right corner, all the diagrams 
contributing to T{c, d,T — 0) are obtained by attaching an external line to a 3-point vertex in the 
diagrams contributing to N'^G{c, d). The diagrams at the first few orders are shown for illustration. 

Then one can obtain Bij (c, d, T) for T / by the shifting the three-point vertex as (See 
Fig. 14) 

Bij{c, d, T) = Bijic d, T = 0) _ _ . (C.8) 
Summation of (C.6) and (C.8) yields 



g^g+\gT 



Bi 



2—G{c,d) + 2—G{c,d)JiiJjj 



- {2XgT + X'gT')i2 + JuJjj) • (C.9) 



J g-*g+XgT 

Prom (C.4) and (C.9) we obtain the Schwinger-Dyson equations for the full propagator 



r + 



d 



d2 



2^0(0, d) 
' 2^0(0, d) 



-2{2XgT + X^gT^) , 



J g-^g+XgT 



{2XgT + X^gT'^' 



(C.IO) 



J g-^g+XgT 



The Schwinger-Dyson equation for the tadpole can be obtained in a similar way. Let 
us first decompose the tadpole T as 



r = r(r) + r(c, d, t) 



(C.ll) 
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The first term, which includes only tadpoles, consists of two diagrams shown in Fig. 13, 
which yield 



r(r) 



1 



(C.12) 



Next we consider the diagrams that contribute to T(c, d,T = 0). Note that every diagram 
other than the one at the lowest order can be obtained from 2PI vacuum diagrams by 
attaching an external line to a 3-point vertex as illustrated in Fig. 15. Thus we have 



f(c, d,T = 0) 



-5-^ (2c + d)+ g-^^Gic d) 



fi + r 



^ + r dg 



(CIS) 



where the first term represents the contribution from the lowest order diagram. Then 
T(c, d, T) with T 7^ can be obtained from T{c, d,T = 0) by the same replacement as in 
Fig. 14, which leads us to 



fic,d,T) 



1 



/U + r 



-g{2c + d)+g—G{c,d) 



(C.14) 



g^g+\gT 



Prom (C.ll), (C.12) and (C.14), we obtain the Schwinger-Dyson equation for the tadpole 



jJi + r 



-g(2c + d)+g^Gic,d) 



.g+XgT f^ + r 



{3XgT^ + X^gT^ 



(C.15) 



Since Eqs. (C.IO) and (C.15) are closed with respect to c, d, T, we can solve them 
order by order once the sum of 2PI diagrams G{c, d) is obtained. 

C.3 Explicit form of the free energy 

Prom the full propagator and the tadpole, one can readily obtain the expectation value of 
the quadratic terms in (C.l) as 



M = ^-^tr02^ = c + Ar^ Mj = (^^tr{J(j)J(l))^=d + XT^. 
Since the free energy ^(//, r) satisfies the differential equations 



dfj,. 



d 



■Hl^,r) =M-Mj , 



(C.16) 



(C.17) 



where fi± = n i r, we can obtain the A-expansion of the free energy by term- by-term 
integration. The integration constant in ^(//, r), which docs not depend on appears 
only at the order of 0(A°), and it can be calculated directly as ^(In AT — In 2). 

The number of 2PI planar vacuum diagrams up to the 3rd order is 13. By evaluating 
them explicitly and by following the above procedure we obtain the explicit form of the 
free energy J^(/x, r) up to the 3rd order as follows. 



= ^ A"J'„(/i,r), 



(C.18) 



n=0 

J^o{^l,r) = ^[21nAr + ln(//+//_) -21n2 
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/ 6 3 1_\ 9(J_ 1 1 \ 

3^ / 108 _ 54 _ 81 _ 17 _ 9 \ 

gf /27 99 47 25 9 \ 

^ 4/i>- 4/x3.;u2 + 4/x2.;u?_ ^ 4/x+/xiy 

y / 9 1_ _ 11 1 9_\ 

^4\l6fiX ^+^' 

56 / 3402 1944 810 351 135 45 9 \ 

y^"^ /1377 1188 1215 306 135 54 27 \ 

g'^f' / 243 621 477 153 81 81 27 \ 

g^ f 9 9 27 9 27 9 9 \ 
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